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Hadronic Interactions at Cosmic Ray Energies 
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General physics of very high energy hadronic interactions is discussed. Special attention is payed to the 
contribution of semihard processes to the interaction dynamics and to the role of parton shadowing and parton 
density saturation. In particular, the implementation of non-linear interaction effects in the QGSJET-II model 
is discussed in detail. The predictions of the model are compared to selected accelerator data, including ones of 
the RHIC collider, and the relation to the calculated extensive air shower characteristics is discussed. Finally, the 
potential of accelerator and cosmic ray experiments for constraining model predictions is analyzed. 



1. INTRODUCTION 

During the last decade a significant progress 
has been achieved in experimental studies of 
high energy cosmic rays (CR) with extensive air 
shower (EAS) techniques, both concerning the 
measurements of the primary cosmic ray (PGR) 
energy spectrum and in determining the compo- 
sition of cosmic radiation in the region of the 
spectral "knee" and at higher energies. This re- 
sulted from further development of the measure- 
ment and data analysis techniques, the latter in- 
cluding thorough EAS simulation studies with 
contemporary Monte Garlo (MG) tools. In EAS 
calculations a special role is played by hadronic 
MG generators which are used for the descrip- 
tion of hadron-air and nucleus-air interactions in 
air showers, performing an extrapolation of cur- 
rent theoretical and experimental knowledge to- 
wards the highest GR energies. In particular, 
the QGSJET model [IE], being based on the 
Gribov's effective Reggeon Field Theory (RFT) 
[3l4l5l6j and the Pomeron phenomenology [718] . 
proved to be very successful in describing air 
shower data obtained by various experimental in- 
stallations. The original versions of this MG gen- 
erator have been developed as MG realizations 
of the Quark-Gluon String (QGS) model [7], in- 
cluding a generalization of the QGS model ap- 
proach for the treatment of nucleus-nucleus col- 
lisions and a description of the fragmentation 



of nuclear spectator part jglTO]. In QGSJET 
this has been supplemented by a phenomenolog- 
ical treatment of semihard processes, which re- 
sult in the production of observed hadron jets of 
comparatively high transverse momenta. Finally, 
in the QGSJET-II model [11] a treatment of 
non-linear interaction effects has been developed, 
based on all-order re-summation of so-called en- 
hanced (Pomeron-Pomeron interaction) RFT dia- 
grams |12I13| . These approaches will be discussed 
below, in comparison with other models and with 
selected accelerator data. 

2. MODEL FRAMEWORK 

2.1. Basic physics 

The general picture for high energy hadronic 
collision is the one of a multiple scattering pro- 
cess, being mediated by multiple parton (quark 
and gluon) cascades proceeding between the two 
hadrons. In the center of mass frame such par- 
ton cascades develop on a much larger time scale 
than the one of a parton-parton scattering, so 
that in the moment of the collision both projec- 
tile and target hadrons are represented by their 
parton clouds. Binary re-scatterings of some of 
these partons violate the coherence of the corre- 
sponding parton chains (consisting of the parton 
"parents" and "pre-parents" ) , which then frag- 
ment into secondary hadrons. Alternatively, the 
coherence may be preserved for some parton re- 
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Figure 1. Proton profile as viewed in soft (a), 
hard (b), semihard (c), and general (d) interac- 
tions at very high energies. 

scatterings, in which case all partons from the 
corresponding chains recombine back to their par- 
ent hadrons without a production of secondaries, 
which gives rise to elastic re-scattering processes. 
A general inelastic collision involves a number of 
elementary hadron production contributions as 
well as multiple elastic re-scatterings. In turn, 
elastic scattering is obtained when only elastic 
sub-processes occur. 

With the energy increasing, the number of el- 
ementary re-scattering processes grows rapidly, 
due to the larger phase space for parton emis- 
sions. In addition, one expects a qualitative 
change in the structure of the underlying par- 
ton cascades. Indeed, at comparatively low en- 
ergies all the partons are characterized by small 
transverse momenta; high pt emissions are sup- 
pressed by the smallness of the corresponding 
running coupling, as{pt). By the uncertainty 
principle, each emission is characterized by a large 
displacement of the produced parton in the trans- 
verse plane, Ab^ ~ l/p?- Thus, with the en- 
ergy increasing further, such "soft" parton cas- 
cades rapidly expand towards larger impact pa- 
rameters, while the density of partons per unit 
transverse area remains small, the hadron look- 
ing "grey", as depicted in Fig. [T] (a). However, 
at sufficiently high energies an important con- 
tribution comes from so-called "semi-hard" and 
"hard" parton cascades, in which some or all 
partons have comparatively high transverse mo- 
menta IH]. There, the smallness of the strong 
coupling as{pt) is compensated by a high par- 
ton density and by large logarithmic ratios of 
the longitudinal and transverse momenta for suc- 
cessive parton emissions. Purely hard cascades, 
which start, e.g., from valence quarks and contain 
only high partons, do not expand transversely. 
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Ab^ ^ \/p1 being small, and lead to an increase 
of parton density in small areas ("hot spots") in 
the transverse plane - see Fig. [l](b), while giving 
a negligible contribution to the total cross section. 
Contrary to that, typical semihard re-scatterings 
are two-step processes: first, parton branchings 
proceed with a small momentum transfer and the 
cascade develops towards larger impact parame- 
ters; next, high pt parton emissions become effec- 
tive, leading to a rapid rise of the parton density 
at a given point in the transverse plane. As a 
result, the region of high parton density extends 
to large impact parameters (Fig. [1] (c)) and the 
contribution dominates in the very high energy 
limit. General hadronic interactions include all 
the mentioned mechanisms; hadrons in high en- 
ergy collisions look as shown in Fig. [1] (d) : there is 
an extended "black" region of high density, dom- 
inated by the semihard processes, and around it 
there is a "grey" region of low density, formed by 
purely soft parton cascading [15]. In the "black" 
region one expects strong non-linear parton ef- 
fects to emerge, which result in the saturation 
of parton densities and in the suppression of soft 
parton emissions [2]. On the other hand, such 
effects are negligible in the "dilute" peripheral re- 
gion. 

How to estimate the relative importance of 
the two regimes? Small peripheral contribution 
would correspond to the "black disc" limit for 
hadron-hadron scattering, with the ratio of elas- 
tic to total cross sections approaching 1/2. The 
relative smallness of the observed (j'^p/cr^pp ratio 
indicates that the "black" central region and the 
peripheral one are yet of comparable sizes [15j . 

In reality, the discussed separation of central 
and peripheral collisions is rather crude, as the 
average parton densities rise gradually with de- 
creasing impact parameter. Thus, there exists 
an important "transition" region of moderately 
large impact parameters, characterized by large 
but not yet saturated parton densities, where the 
contributions of both soft and semihard processes 
are of equal importance, and where non-linear 
parton effects provide sizable corrections. In fact, 
it is this transition region which is expected to 
give the dominant contribution to the interaction 
characteristics relevant for FAS physics. 
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Figure 2. A general contribution to hadron- 
hadron scattering amplitude. Elementary scat- 
tering processes (vertical thick lines) are de- 
scribed as Pomeron exchanges. 

2.2. Pomeron formalism 

As discussed above, general hadronic collisions 
necessarily involve emissions of soft low pt par- 
tons, which prevents one from applying the per- 
turbative QCD formalism. However, for a num- 
ber of key characteristics, like total and elastic 
cross sections or probabilities of various "macro- 
scopic" configurations of inelastic interactions, 
the knowledge of the microscopic parton picture 
is not of extreme importance; the RFT allows 
one to calculate the quantities of interest based 
on the knowledge of the elastic amplitude for an 
"elementary" scattering process. For example, 
in the QGS model [7j hadron-hadron scattering 
is described as a multiple exchange of composite 
objects - Pomerons, as shown in Fig. [21 Each 
Pomeron represents a microscopic parton cas- 
cade whose precise description is not necessary at 
this stage. Rather one employs the theoretically- 
motivated ansatz for the corresponding scattering 
amplitude: 



Qp(0)-1 



i^Ri+Ri+c'm Ins) 



ai,(0)lns 



-,(1) 



which is characterized by a power-like energy rise 
and by a logarithmically increasing slope; the pa- 
rameters ap(0), ap(0) are the intercept and the 
slope of the Pomeron Regge trajectory and 7a, 
- the residue and the slope for the Pomeron- 
hadron a vertex [7)B\ . 

It is important to take into account the con- 
tributions of both elastic and inelastic interme- 
diate hadron states between Pomeron emissions, 
such that the Pomeron-hadron vertex becomes a 
matrix corresponding to the transitions between 
those states. To diagonalize this matrix one con- 



Figure 3. Typical inelastic interaction contains a 
number of elementary production processes, de- 
scribed by cut Pomerons (thick broken lines in the 
Figure), and any number of elastic re-scatterings 
- uncut Pomeron exchanges. 

siders any hadron to be represented by a super- 
position of a number of elastic scattering eigen- 

states, I a) = J2j \J ^a'' \ja), which are charac- 
terized by different absorption in the scattering 
process, — Xa^ [6J. Here Ca^ and Aq^ 
are the weights of the eigenstates and the rela- 
tive strengths for their coupling to the Pomeron; 



EjCi''^ = 1- EjCk'^X'a" = 1. Then, the to- 
tal cross section is obtained as the sum of partial 
scattering contribitions, being expressed via the 
so-called Pomeron eikonal X^A^^ ^) ~ ^^fadi^' ^) 

'^TXs) = 2/^^55: C(^-)cf(l-e^^i^'^i^'>^-(-''))(2) 

To obtain cross sections for various final states, 
one makes use of the optical theorem, which re- 
lates the total sum of contributions of all final 
states to the imaginary part of the elastic ampli- 
tude for hadron-hadron scattering, hence, to the 
contributions of various unitarity cuts of elastic 
scattering diagrams. The so-called Abramovskii- 
Gribov-Kancheli cutting rules [12] state that only 
certain classes of cut diagrams are important in 
the high energy limit and allow one to relate such 
contributions to particular final states of interest. 
For example, the cross sections for having n si- 
multaneous elementary production processes are 
described by the diagrams with n "cut" Pomerons 
and any number of uncut ones, as shown in Fig.[3l 
whereas elastic and diffractive cross sections are 
obtained cutting the diagrams of Fig. [5] between 
the Pomerons, with no one being cut, and se- 
lecting elastic or diffractive intermediate hadron 
states in the cut plane. Thus, the cross sections 
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for the n cut Pomerons process and for the pro- 
jectile diffraction dissociation read 




From (j4|) foUows that a higher diffraction cross 
section is obtained for more asymmetric couphngs 
Xa(j) of different eigenstates to the Pomeron. In 
turn, this resuhs in a smaller total cross section 
and larger fluctuations in the number n of elemen- 
tary production processes (cut Pomerons), as one 
can see from ([5H31). 

Having obtained partial probabilities of vari- 
ous inelastic final states, one has yet to describe 
particle production for each elementary inelastic 
process, which leads us back to the physics of the 
underlying parton cascade. However, one can em- 
ploy here the string picture of the hadronization, 
based on the color structure of the correspond- 
ing final states. Namely, one assumes that each 
cut Pomeron process induces a color exchange 
between parton constituents of the interacting 
hadrons, such that color strings are stretched be- 
tween them. With the two hadrons flying apart, 
these strings break up and hadronize, which is 
described by string fragmentation procedures, ft 
is noteworthy, that in the QGS model the pa- 
rameters for constituent parton momentum dis- 
tributions and the ones for string fragmentation 
are not adjustable ones, being expressed via the 
characteristics of known Regge trajectories |7|17j . 

Remarkably, the described scheme can be gen- 
eralized to hadron-nucleus and nucleus- nucleus 
interactions in a parameter-free way, both for 
cross section calculations and concerning the de- 
scription of particle production |10| . The only 
new input are nuclear density profiles, with the 
corresponding parameters being fixed by nuclear 
form factor measurements |18] . individually for 
each nuclear type, as discussed in ilOj. 




Figure 4. A general Pomeron is the sum of the 
"soft" and the "semihard" ones - correspondingly 
the 1st and the 2nd graphs in the r.h.s. 

2.3. Semihard processes 

As discussed in Section [^TTl at high enough en- 
ergies a significant contribution to hadron-hadron 
scattering comes from semihard parton processes, 
where a part of the underlying parton cascade 
develops in the high pt region. To provide a mi- 
croscopic treatment of the corresponding physics, 
one can apply the phenomenological Pomeron de- 
scription to the low virtuality — Pt < Qq) 
part of the parton cascade, while treating the high 
Pt parton evolution using pQCD techniques, Qq 
being some chosen virtuality cutoff for QCD be- 
ing applicable. Then, an elementary scattering 
process is described as an exchange of a "general 
Pomeron" , which is a sum of the "soft" and the 
"semihard" ones, the latter being represented by 
a piece of QCD ladder sandwiched between two 
"soft" Pomerons [Il2l20l21j . as shown in Fig. H 
As discussed in |20|21) . the upper and the lower 
half of the semihard Pomeron define parton (sea 
quark and gluon) momentum and impact pa- 
rameter distributions in the projectile and target 
hadrons; the parameters for the coupling between 
the soft Pomeron and the ladder can be fixed 
by the data on hadron structure functions (SFs). 
The approach allows one to treat the semihard 
processes within the Pomeron scheme described 
in Section 12. 2| the interaction cross sections are 
defined by the usual formulas, like ([2H11), with the 
Pomeron eikonal being a sum of the soft and the 
semihard ones, xldi^' ^) = Xad'^ (s, b)+x'ad{s, b). 
In turn, particle production procedure includes 
an explicit treatment of the high pt parton cas- 
cade, with the strings being formed between the 
produced final partons [2l20l21j . 
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2.4. Non-linear effects 

Developing a model for high energy hadronic 
interactions, one inevitably faces the problem of 
treating non-linear effects, connected to parton 
shadowing and saturation. Indeed, describing 
the interaction as a superposition of a number of 
re-scattering processes, mediated by parton cas- 
cades, one has to consider the case when such cas- 
cades overlap in the corresponding phase space 
and influence each other. Such effects are ex- 
pected to be extremely important at very high 
energies and small impact parameters, i.e. in the 
"black" region of high parton densities, where 
they lead to the parton density saturation [M] 
and to significant reduction of secondary particle 
production. However, non-linear effects start to 
be efficient already at comparatively low energies 
and large impact parameters, the experimental 
indication being the rapid energy rise of the high 
mass diffraction cross section in the ISR energy 
range [22], the latter being just one of a number 
of manifestations of non-linear parton dynamics. 
Therefore, one has to develop a coherent descrip- 
tion of the corresponding physics over a wide dy- 
namic range, rather than restrict himself with a 
treatment of the saturation region. 

Treating independent parton cascades effec- 
tively as Pomeron exchanges, the corresponding 
non-linear effects are described in the RFT as 
Pomeron-Pomeron interactions |23|24|25] . There, 
the main technical difficulty is to perform a re- 
summation of contributions of the underlying en- 
hanced Pomeron graphs, as more and more dia- 
grams of compicated topologies come into play 
at higher energies. A re-summation method 
has been worked out recently [12ll3j and imple- 
mented in the QGSJET-II model [11] ■ The basic 
assumption of the approach was that Pomeron- 
Pomeron coupling proceeds via parton processes 
at comparatively small virtualities < Qo ^^d 
can be described using phenomenological multi- 
Pomeron vertices of eikonal type PJJ. A rea- 
sonable consistency with relevant experimental 
data has been obtained using a fixed energy- 
independent Qo- cutoff0 neglecting parton shad- 

'^The parameter set in 1131 is different from the default 
QGSJET-II settings. While a general agreement with 
measured proton-proton cross section and proton SFs can 
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pseudorapidity eta 

Figure 5. Pseudorapidity distributions of charged 
secondary particles in Au-Au collisions of differ- 
ent "centralities" ; QGSJET-II results (lines) are 
compared to BRAHMS data [27] . 

owing effects for > Qq. In particular, infer- 
ring the basic parameter of the scheme, the triple- 
Pomeron coupling, from HERA hard diffraction 
data, one obtains a satisfactory agreement with 
observed hadronic diffraction [l3] . 

To describe particle production, one has to con- 
sider unitarity cuts of enhanced Pomeron dia- 
grams and, for any configuration of cut Pomerons, 
to perform a full re-summation of uncut ones 
|13I26| . The big advantage of the developed pro- 
cedure is that the solutions are obtained in the 
form of recursive equations, which can be imple- 
mented in a MC model and allow one to generate 
various configurations of interactions, including 
diffractive ones, in an iterative fashion [26] . 

The formalism seems to be adequate for the 
description of hadronic collisions in the periph- 
eral and "transition" regimes. However, in the 
"black" region of high parton densities one may 
expect an important contribution of the "hard" 
(k^l > Qo) Pomeron-Pomeron coupling, which 
would provide additional screening corrections. A 
reasonable agreement of QGSJET-II with RHIC 
data on central heavy ion collisions (see Fig. [5]) 
indicates the smallness of such effects. However, 
the situation may change at much higher energies. 

be obtained for Qq as small as 1 GeV^, the description 
of particle production requires a higher cutoff, the default 
QGSJET-II value being 2.5 GeV^. This may be consid- 
ered as an indication for significant parton shadowing in 
the virtuality range 1 2.5 GeV^. 
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2.5. Nuclear fragmentation 

Performing simulations of nucleus-induced air 
showers, one has to pay attention to the frag- 
mentation of nuclear spectator part, as it has a 
significant influence on the predicted EAS fluctu- 
ations [To]. Nuclear fragmentation in QGSJET 
is strongly correlated with the interaction treat- 
ment; it is assumed that the excitation of the nu- 
clear spectator part is proportional to the num- 
ber of inelastically "wounded" nucleons. In case 
of comparatively low excitations, corresponding 
to peripheral interactions, the de-excitation pro- 
ceeds via an evaporation of a certain number of 
nucleons or alpha-particles. Highly excited nu- 
clei undergo multi-fragmentation, which is de- 
scribed as a percolation process: final fragments 
are formed by the clusters of nearby nucleons [lOj . 
The scheme provides a reasonably good descrip- 
tion of the mass yield dependence for all fragment 
masses, including intermediate ones, which is typ- 
ically not the case if evaporation mechanism only 
is included (e.g., in [H])- Importantly, the rele- 
vant model parameters have been fixed with the 
data on nuclear fragmentation at energies above 
2 AGeV, where the measured fragment yields are 
already stable and do not have any significant en- 
ergy dependence. Correspondingly, the treatment 
can be safely extrapolated to cosmic ray energies. 

3. DISCUSSION OF MODEL RESULTS 

3.1. Cross sections and inelasticity 

The longitudinal air shower development de- 
pends on a number of characteristics of hadron- 
air interactions. While the largest effect comes 
from the inelastic proton-air (nucleus-air) cross 
section, which defines the position of the first in- 
teraction of the primary particle, pion-air and 
kaon-air cross sections and the inelasticities of 
hadronic interactions influence the shape of the 
characteristic EAS profile. In particular, cross 
sections of diffraction dissociation processes play 
here a special role. The diffraction of the tar- 
get nucleus results in a negligible energy loss of 
the initial particle, secondary hadrons being pro- 
duced in the target fragmentation region, and is 
essentially equivalent to a reduction of the inelas- 
tic cross section. The projectile hadron diffrac- 



Table 1 

Inelastic and non-diffractive hadron-carbon cross 
sections (in mb) at 200 GeV lab. energy. 



QGSJET-II 238 224 167 

Experiment 237 ± 2 [22] 225 ± 7 [3D] 171 ± 5 [50] 



tion, being predominantly a coherent process, the 
target nucleus remaining intact with a high prob- 
ability [TU], results in a small energy loss of the 
leading particle, thus sizably influencing the over- 
all inelasticity. It is noteworthy, that all the 
above-mentioned quantities mainly depend on the 
model description of hadronic collisions in the 
peripheral and "transition" regions. Indeed, at 
small impact parameters, hadrons already look 
"black", see Fig. [T] (d), with the interaction pro- 
file cr™f^j.(s, b) (the probability of inelastic inter- 
action at some impact parameter b for a given 
energy s) being very close to unity, which still 
remains so when non-linear corrections to par- 
ton dynamics, e.g., parton saturation effects, are 
taken into account. This region does not con- 
tribute substantially to diffraction cross sections, 
as any rapidity gap, produced in some elementary 
production process, is covered by secondary par- 
ticles created in other production sub-processes. 
Finally, the inelasticity can only weakly depend 
on the central region treatment. Central colli- 
sions involve a large number of production pro- 
cesses and lead in all cases to large energy losses of 
leading particles. Contrary to that, modifications 
of the treatment of the peripheral and the "tran- 
sition" regions make a sizable effect on the pre- 
dicted a"f|jij.(s, 5) behavior, thus influencing the 
obtained values of cr™l\i^ and o-^'^ain and change 
the average number and fluctuations of the num- 
ber of production sub-processes, which leads to 
large variations of the inelasticity. 

As already mentioned, the generalization of 
the formalism, discussed in Sections I2.2H2.41 for 
hadron-nucleus and nucleus-nucleus interactions 
proceeds in a parameter-free way, which applies 
also to cross section calculations. For the lat- 
ter, one obtains an excellent agreement with ex- 
perimental data, some values given in Table [1] 
The predictions of the SIBYLL 2.1 model |31I32] 
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for a'^^l-^. are close to the ones of QGSJET-II 
in the coUider range but rise faster with energy. 
In principle, this may come from the fact that 
non-linear effects are introduced in SIBYLL for 
semihard processes only, being neglected for the 
"soft" component. As the latter plays the cru- 
cial role in peripheral collisions, which largely de- 
fine the cross section behavior, one may expect 
that SIBYLL overestimates cr^fair ^^'^ very 
high energy range. Nevertheless, model predic- 
tions depend also on the assumptions concern- 
ing hadron form factors (parton distributions in 
the transverse plane) and on the calibration to 
proton- proton data, as discussed in [33 . Present 
uncertainties concerning the energy dependence 
of cr^"f!^jj. will be drastically reduced after the mea- 
surements of (Tp°* at the LHC. 

Model results for the inelasticity ^i^/j'lfair '^^ 
hadron-air interactions are only constrained by 
fixed target accelerator measurements. Practi- 
cally in all models the inelasticity rises with en- 
ergy, due to the increase of the number of el- 
ementary production processes. However, the 
spread in the predicted i^^^faij. is rather large, 
around 20% in the high energy range. Additional 
constraints may come from the measurements of 
leading neutron spectra in proton-proton interac- 
tions by the LHCf experiment 34J. 

3.2. Muon component of air showers 

An important topic are model predictions for 
the CR muon component. Here one has to distin- 
guish between the results for inclusive muon spec- 
tra and for the muon content of air showers of a 
given primary energy. The former are dominated 
by single interactions of primary protons of ener- 
gies in average only an order of magnitude higher 
than the ones of the measured muons. Due to the 
steepness of the primary CR spectrum, the cor- 
responding results are very sensitive to the shape 
of the forward pion and kaon spectra in proton- 
air collisions. The characteristic quantities are 
the so-called Z- factors some values given 

in Table [2] for QGSJET-II and SIBYLL 2.1. 
Clearly, the two models predict different energy 
dependences for : one observes a rather pre- 
cise Feynman scaling in SIBYLL, which is sup- 
ported by inclusive muon flux measurements |36| , 



Table 2 



Pion and kaon Z-factors as predicted by the 
QGSJET-II and SIBYLL 2.1 models. 



Eo, TeV 


model 






Zk+ 


0.1 


QGSJET-II 
SIBYLL 2.1 


0.043 
0.036 


0.035 
0.026 


0.0036 
0.0134 


10 


QGSJET-II 
SIBYLL 2.1 


0.033 
0.037 


0.028 
0.029 


0.0034 
0.0097 



and a noticeable scaling violation in QGSJET- 
II. Another difference is the higher values of 
in case of SIBYLL, which is due to the harder 
kaon spectra in the model. 

The EAS muon content is formed during a 
multi-step hadronic cascade process and mainly 
depends on the total multiplicity of hadron-air 
collisions, the shape of the forward pion spectra 
being of secondary importance. Model calibra- 
tion is mainly performed at fixed target energies; 
in particular, both QGSJET-II and SIBYLL 
2.1 appear to be consistent with recent data of 
the NA49 Collaboration [37j . Concerning the pre- 
dicted multiplicity, models agree with each other 
in the collider range; at the highest CR ener- 
gies the characteristic differences reach a factor 
of three. Most likely, the situation will not im- 
prove significantly with the start of the LHC, as 
the spread in model results for the multiplicity of 
proton-proton interactions is yet at 10% level in 
the corresponding energy range. 

However, one should not expect large differ- 
ences between model predictions for the EAS 
muon number N^. In particular, QGSJET-II 
and SIBYLL 2.1 results for differ by only 
10 15% in the high energy range. As shown in 
[38] , changing the multiplicity of proton-air col- 
lisions by a factor of two, one obtains less than 
10% modification of the predicted N^j,. Of course, 
pion- air multiplicity is of greater importance here. 
However, in QGSJET-II one has little freedom 
to modify the latter. The only new parame- 
ters, defining the transition from proton-proton 
to pion-proton case, are the residue and the 
slope of the Pomeron-pion coupling, which are 
reliably fixed by data on pion-proton total cross 
section and elastic scattering slope. The low x 
behavior of the gluon and sea quark momentum 
distributions in pion are described by the same 
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soft Pomeron asymptotics as in the proton case 
and the corresponding normahzation is defined by 
the momentum sum rule, the valence quark distri- 
butions being fixed by measurements. It is there- 
fore surprising that the new EPOS model [39] 
predicts significantly higher iV^, which is almost 
a factor of two in access of the QGSJET-II re- 
sults at the highest CR energies [40] . Most likely, 
this is connected to the treatment of pion-air in- 
teractions in EPOS. The picture can be validated 
by the KASCADE experiment, where the muon 
component is measured directly, the correspond- 
ing analysis being sensitive to 10% variations in 
the predicted muon number |41| . 

The author acknowledges fruitful discussions with 

H. J. Drescher, R. Engel, and M. Strikman. 
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